Semiclassical approach to the ac-conductance of chaotic cavities 
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We address frequency-dependent quantum transport through mesoscopic conductors in the semi- 
classical limit. By generalizing the trajectory-based semiclassical theory of dc quantum transport to 
the ac case, we derive the average screened conductance as well as ac weak-localization corrections 
for chaotic conductors. Thereby we confirm respective random matrix results and generalize them 
by accounting for Ehrenfest time effects. We consider the case of a cavity connected through many 
leads to a macroscopic circuit which contains ac-sources. In addition to the reservoir the cavity 
itself is capacitively coupled to a gate. By incorporating tunnel barriers between cavity and leads 
we obtain results for arbitrary tunnel rates. Finally, based on our findings we investigate the effect of 
dephasing on the charge relaxation resistance of a mesoscopic capacitor in the linear low-frequency 
regime. 
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I. INTRODUCTION 

In contrast to dc-transport experiments, the applied 
external frequency w of an ac-driven mesoscopic struc- 
ture provides a new energy scale hw that permits one to 
access further properties of these systems, including their 
intrinsic charge distribution and dynamics. 

The interest in the ac-reponse of mesoscopic conduc- 
tors goes back to the work of Pieper and Priced on the 
dynamic conductance of a mesoscopic Aharonov-Bohm 
ring. This pioneering work was followed by several exper- 
iments ranging from photon-assisted transport to quan- 
tum shot noise^ii^'^'^'i. More recently, the ac-regime has 
been experimentally reinvestigated achieving the mea- 
surement of the in and out of phase parts of the ac- 
conductance^ and the realization of a high-frequency sin- 
gle electron source^. Moreover, the recent rise of interest 
in the full counting statistics of charge transfer has led 
to a reexamination of the frequency noise spectr a^°i^^i^^ . 
This experimental progress has since triggered renewed 
theoretical interest in time dependent mesoscopic trans- 

p □1- ^13.14.15.16.17 

One way to tackle the ac-transport problem is to start 
from linear response theory for a given potential distri- 
bution of the sampl o^^i^^i^° . This involves the difficulty 
that, in principle, the potential distribution and more 
precisely its link to the screening is unknown. Another 
approach consists of deriving the ac-response to an ex- 
ternal perturbation that only enters into quantities de- 
scribing the reservoirs. Such approachs were initiated 
by Pastawski^i within a non-equihbruium Green func- 
tion based generalized Landauer-Biittiker formalism, and 
then the scattering matrix formalism of a time-dependent 
system was developed by Biittiker et al.— >2^. Since the 
energy is in general no longer conserved for an ac-bias, 
the formalism is based on the concept of a scattering ma- 
trix that depends on two energy arguments^ or equiv- 
alently on two times^. Fortunately, when the inverse 
frequency is small compared to the time to escape the 



cavity, the ac-transport can be expressed in terms of the 
derivative of the scattering matrix with respect to en- 
ergy2£. In this article we start from the time dependent 
scattering matrix formahsm and Hmit our investigations 
to open, classically chaotic ballistic conductors in the low- 
frequency regime^l. 

For ac-transport we calculate the average correlator of 
scattering matrices S(i?) at different energies E. For this 
we need to know the joint distribution of the matrix el- 
ements Sai3;ij at different values of the energy or other 
parameters. (We label the reservoirs connected to the 
conductor by a greek index and the mode number by a 
latin index.) To our knowledge a general solution to this 
problem does not yet exist for chaotic systems. How- 
ever, in the Hmit of a large number of channels, the first 
moments of the distribution Sa/3;ij{E)Sl^p.-j{E') were de- 
rived using both semiclassical method a^^i^^ and various 
random matrix theory (RMT) based method92^>^>^>2^. 
Although the ac-transport properties of ballistic chaotic 
systems seem to be well described by the RMT of trans- 
porlj^, we develop a semiclassical approach for three rea- 
sons: First, this allows us to confirm the random matrix 
prediction by using a complementary trajectory-based 
semiclassical method. Second, the energy dependence 
in the random matrix formalism was introduced by re- 
sorting to artificial models such as the "stub model"—. 
While being powerful, this treatment is far from micro- 
scopic or natural. The third and strongest reason is to go 
beyond the RMT treatment and investigate the crossover 
to the classical limit. Similarly as for the static case 
RMT is not applicable in this regime. As first noticed by 
Aleiner and Larkin^, ballistic transport is characterized 
by a new time scale, known as the Ehrenfest time tet^h^, 
that controls the appearance of interference effects. The 
Ehrenfest time corresponds to the time during which a 
localized wavepacket spreads to a classical length scale. 
Typically, in open chaotic systems two such lengths are 
relevant, the system size L and the lead width W. We 
can thus define an Ehrenfest time associated with each 



2 



one^i^, the closed-cavity Ehrenfest time, 

r^' = A-Mn[L/AF], (1) 
and the open-cavity Ehrenfest time, 

r°P = A-Mn[M^VAFi], (2) 

where A is the classical Lyapunov exponent of the cavity. 

Although the success of the semiclassical method 
(beyond the so-called diagonal approximation, see be- 
low) to describe quantitatively universal and non 
universal dc-transport properties is now clearly es- 
tablished2S>29dodid2d3d4d5d6dld8d9, the correspond- 
ing semiclassical understanding of frequency dependent 
transport is far less developed. Based on an earlier 
semiclassical evaluation of matrix element sum rules by 
Wilkinson^ and a semiclassical theory of linear response 
functions^, a semiclassical approach to the frequency- 
dependent conductivity within the Kubo-formalism led 
to an expression of the ac-(magneto-) conductivity o'(lj) 
in terms of a trace formula for classical periodic orbits^. 
Closely related to this evaluation of a{uj) is the problem 
of frequency-dependent (infrared-) absorption in ballistic 
mesoscopic cavities which has been treated semiclassi- 
cally in Ref. [5l|. Peaks in the absorption could be as- 
signed to resonance effects when the external frequency 
iv corresponds to the inverse periods of fundamental pe- 
riodic orbits in the cavity. Ref. [s^ contains a first, 
tr-model based approach to weak localization effects in 
the ac-Kubo conductivity, where the findings were inter- 
preted in a quasiclassical trajectory picture (beyond the 
diagonal approximation). We note also that the semi- 
classical treatment of the product of scattering matrices 
S{E) at different energies, has been investigated in dif- 
ferent context such as the Ericson fluctuations^ and the 
time delay^, however without considering the Ehrenfest 
time dependence. 

The outline of this article is as follows: In Section lUl 
we introduce our model to treat the system of interest 
namely a quantum dot under ac bias, and recall some ba- 
sic results about conservation laws in presence of a time 
dependent field. In Sect. [Till we present the method used 
to treat screening, which is based on a self-consistent ap- 
proach developed by Biittiker et ali^. The admittance, 
i.e. the ac-conductance, is then calculated semiclassically 
for the particular case of strong coupling to the leads 
(transparent contact) in Sect. II VI where we illustrate our 
result by treating the time dependence of a pulsed cavity. 
We generalize the method to cope with arbitrary tunnel 
rates in Sect. |Vl and finally we use our general results 
to investigate dephasing effects on the charge relaxation 
resistance of a mesoscopic capacitor in Sect. I VII 

II. THE MODEL 

We consider a ballistic quantum dot, i.e. a two- 
dimensional chaotic cavity coupled to M electron reser- 
voirs via M leads. Each lead a has a width Wa and 



Wo(aj),2'o(w) 




Figure 1: Two dimensional chaotic cavity with M leads and 
one gate 0. Each lead a has a width Wa and is coupled to 
a reservoir at potential Uaii^) and current 3a{(jj)- Each tun- 
nel barrier is characterized by the set of transmission prob- 
abilities Ta = {ra,i,--- ,ra_]v„}. The gate and the sam- 
ple are capacitively coupled, which leads to a gate current 
3o(w) = -iujC[Uo{u}) - Uiuj)]. 



is coupled to the cavity through a tunnel barrier (see 
Fig. [T|) . In addition to the treatment of Ref. [i^ we as- 
sign a particular tunnel probability to each lead mode. 
The tunnel barrier is thus characterized by a set of trans- 
mission probabilities, Fq, = {r„^i, • • • jFa^Ar^}, with Na 
the maximum mode number of lead a. The chaotic dot is 
additionally capacitively coupled to a gate connected to 
a reservoir at voltage 'Uo(w), from which a current 3o{uj) 
flows. This capacitive coupling with the gate is taken 
into account via a geometrical capacitance C^^i^^i^^ . 

We further require that the size of the contact is much 
smaller than the system size L, but still semiclassically 
large, 1 ^ Na ^ L/Ap. This requirement ensures that 
the particle spend enough time inside the cavity to expe- 
rience the chaotic dynamics. 

As usual for such mesoscopic structures we need to dis- 
tinguish between quantum and classical time scales. On 
the quantum side we have already introduced the Ehren- 
fest times {t^ , r^' ) in Eqs. (|ll2p . while another time scale 
is the Heisenberg time th, the time to resolve the mean 
level spacing of the system. On the classical side the time 
of flight Tf between two consecutive bounces at the sys- 
tem cavity wall is relevant. In most ballistic systems or 
billiards we have T[ ~ A~^. Another relevant time scale 
is the ballistic ergodic time Torg which determines how 
long it takes for an electron to visit most of the available 
phase space. However, as we deal with transport proper- 
ties, a further important time scale is the dwell time td, 
the average time spent in the cavity before reaching the 
contact, we have Tn/Tag 3> 1. The related escape rate 
therefore satisfies 

M 
Q— 1 i—1 

For small openings which we consider here, we have 
Atd > 1. 
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The ac-transport properties of such a mesoscopic sys- 
tem are characterized by the dimensionless admittance 

go,p{uo) = Go,p{uj)IGo = G^^dIo,{uj)/dUp{u), (4) 

with Go = dse^/h, where ds = 1 or 2 in the absence 
or presence of spin degeneracy. In this study we hmit 
ourselves to the coefficients gap{uj) with a,(3 = 1, • • • , M 
where the coefficients denoting the gate are determined 
by current conservation and the freedom to choose the 
zero point of energy^^, 



M 

E 

Q=0 



M 



9ap{i^) = ^gap{i^) = 0. 



(5) 



We note that Eq. Q is a straightforward consequence 
of the underlying gauge invariance. Owing to the con- 
servation of charge, the total electric current fulfills the 
continuity equation 



dt ' 



(6) 



where p is the charge density and jp the particle current 
density. For dc-transport, the charge density is time in- 
dependent and so we have V • jp = 0. Thus the sum of all 
currents that enter into the dot is always zero. Moreover 
the current properties must remain unchanged under a si- 
multaneous global shift of the voltages of the reservoirs. 
These conditions imply the well know unitarity of the 



scattering matrix^. 



(7) 



For ac-transport, the product of scattering matrices 
at different energies no longer obey a similar prop- 



erty^>^>^>^ i.e. 



(8) 



indeed this inequality expresses the fact that, due to the 
possible temporary pile up of charge in the cavity, the 
particle current density no longer satisfies V • jp = 0. 
However one can instead use the Poisson equation 



V • D = p. 



(9) 



where D = —eoVf with the electric potential, to define 
the total electric current density which satisfies V • j — 0, 
as a sum of a particle and a displacement current: 



9D 

'dt' 



(10) 



In order to find j one needs to know the electrical field 
D. In general its calculation is not a trivial task because 
the intrinsic many-body aspect of the problem makes the 
treatment of the Poisson equation ([9]) tricky, especially 



if it is necessary to treat the particle and displacement 
current on the same footing. 

In this work we shall adopt the approach of Ref. jl^l to 
simplify the problem. In this approach the environment 
is reduced to a single gate, the Coulomb interaction is 
described by a geometrical capacitance C, and the two 
currents are treated on different footing; the particle cur- 
rent is calculated quantum mechanically via the scatter- 
ing approach, while the displacement current is treated 
classically via the electrostatic law (Eqs. (|6l9p ). This sim- 
plification will permit us below to re-express the Poisson 
equation ^ to obtain the simplest gauge invariant the- 
ory that takes care of the screening. We emphasize that 
even though our model could be thought of as oversim- 
plified it has the advantage of being able to probe the 
effects due to the long range Coulomb interaction. In- 
deed, for non-interacting particles it is possible to treat 
the dot and the gate via two sets of uncorrelated conti- 
nuity equations. The Coulomb interaction removes this 
possibility, and we need to consider the gate and dot as 
a whole system. 



III. EXPRESSION FOR THE ADMITTANCE 

The method to compute the admittance proceeds in 
two steps^: First the direct response (particle current) 
to the change of the external potential is calculated un- 
der the assumption that the internal potential U{uj) of 
the sample is fixed. This leads to the definition of the 
unscreened admittance g'^^p{uj). Second, a self-consistent 
procedure based on the gauge invariance (current con- 
servation and freedom to choose the zero of voltages) is 
used to obtain the screened admittance gapiuj). 

The unscreened admittance reads^^ 



glpiu)^ JdE 

X Tr Sapla ^ 



'q/3 



(11) 



Q/3 
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where f{E) stands for the Fermi distribution, Sap is the 
Na X Np scattering matrix from lead {3 to lead a, and Iq, 
is an Na x Na identity matrix. Under the assumption 
that U(a;) is spatially uniform, the screened admittance 
gafs{^) is straightforward to obtain^. For sake of com- 
pleteness we present here only the outline of the method 
and refer to Ref. [1^ for more details. 

On the one hand the current reponse at contact a is 



3a (w) = Go 



M 

E 

13=1 



glp{u:)Up{u;) + gULo)U{Lo) 



(12) 



where 5^0 ('^) is the unknown internal reponse of the 
mesoscopic conductor generated by the fiuctuating po- 
tential 11(0;). On the other hand the current induced at 
the gate is 



%{uj) = ^iuC[Uo{uj)-U{uj)]. 



(13) 
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Gauge invariance permits a shift of —IKuj) and provides 
an expression for the unknown internal response, 



M 
/3=1 



(14) 



Then current conservation, X)a=i-'a(^) + "^oi^) ~ 0, 
yields the result of the screened admittances^. 



(15) 



In the self-consistent approach used to obtain Eq. (|T5|) . 
the only electron-electron interaction term that has been 
considered is the capacitive charging energy of the cavity. 
This implies that we should consider a sufficiently large 
quantum doli^. We note that, using a 1/A^-expansion, 
the self-consistent approach above was recently formally 
confirmed in Ref. [59|. Moreover, Eq. lfT5|) can be gener- 
alized to non-equilibrium problems, using Keldysh non- 
equilibrium Green functions^. 

In the next sections we present the semiclassical evalu- 
ation of Eq. (flTI) in the zero temperature limit (including 
finite temperature is straightforward) . For reasons of pre- 
sentation we first give the semiclassical derivation for the 
transparent case in Sect. ITVl and then we explore the 
general case in Sect.O In Sect. I VII we present an appli- 
cation of the screened result for tunnel coupling, when we 
compute the relaxation resistance of a mesoscopic chaotic 
capacitor. 



IV. SEMICLASSICAL THEORY FOR THE 
ADMITTANCE 



paths 7 connect Xq = {xo,Pxq) (on a cross section of 
lead P) to X = {x,px) (on a cross section of lead a). 
Each path gives a contribution oscillating with action 
(including Maslov indices) evaluated at the energy 
zL hw/2 and weighted by the the complex amplitude 

. This reduces to the square root of an inverse element 
of the stability matrix^, i.e. = \ {dpj.o/dx)^\^ . 

We insert Eq. lfT7|) into Eq. lfT6|) and obtain double 
sums over paths 7, 7' and lead modes \i), \j). The sum 
over the channel indices is then performed with the semi- 
classical approximation^, J2i=i{^o\'i'){i\xQ) Ri 6{x'q — xo), 
and yields 

9a0i^) ~ NaSa/3 = - dxo /da: ^ 



A A* 

2iih 



7,7 



Here, 



(18) 



5S{Ef,uj) = S^{xo,x]Ef + — ) - Sy{xi),x]EF ^). 

(19) 

As we are interested in the limit hw <C Ey , we can expand 
5S{Ey,oj) around E-p. The dimensionless ac-conductance 
is then given by 



dxo / dx 2^ 

3 



2TTh 



(20) 



X exp 



L5S{EY)+'-^{t,+t,,) 



where 5S{Ey) = S~f{xQ,x; E-p) — S^'{xq,x] Ep) and 
(t~^i ) is the total duration of the path 7 (7'). Eq. (|20|) is 
the starting point of our further investigations. 



Drude Admittance 



A. Semiclassical approximation 

We first consider the multi-terminal case assuming 
transparent barriers, i.e. Ta,i = 1, V(a,i). In the limit 
k^T — > the unscreened admittance, Eq. (fTTj) . reduces 
to 



S.,(i^F + ^)Sl,(i?F-^) 



(16) 

Semiclassically, the matrix elements for scattering pro- 
cesses from mode i in lead {3 to mode j in lead a rea d^^i^^ 



Saf}:,ji{Ep ± — ) 



(17) 



,{i\x){xo\i) 



where \i) is the transverse wave function of the i-th mode. 
Here the (or x) integral is over the cross section of the 
/?th (or ath) lead. At this point Sq^j is given by a sum 
over classical trajectories, labelled by 7. The classical 



We are interested in quantities arising from averaging 
over variations in the energy or cavity shapes. For most 
sets of paths, the phase given by the linearized action dif- 
ference 5S{Ey) will oscillate widely with these variations, 
so their contributions will average out. In the semiclas- 
sical Hmit, the dominant contribution to Eq. l(20l) is the 
diagonal one, 7 
and gives 



7', which leads to = t^i , SS{Ep) ~ 



glf{u;)=N^6^^^ JdxoJdxY.^-^e^'^'^- (21) 

In the following we proceed along the lines of Ref. [43 |. 
The key point is the replacement of the semiclassical am- 
plitudes by their corresponding classical probabilities. To 
this end we use a classical sum rule valid under ergodic 



assumptiona*^. 



(22) 



/ dt d9od9e'^'pFCOs{0o)P(X,Xo;t)[---U„. 

Jo J~tt/2 
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Encounter 




Figure 2: A semiclassical contribution to weak localization 
for a system with strong (transparent) coupling to the leads. 
The two paths follow each other closely everywhere except 
at the encounter, where one path (dashed line) crosses itself 
at an angle e, while the other one (full line) does not (going 
the opposite way around the loop). The cross-hatched area 
denotes the region where two segments of the solid paths are 
paired (within Wa ^ ~ of each other) 



In Eq. (|22|) . cos(0o) is the initial momentum along the 
injection lead and P(X,Xo;i) the classical probability 
density to go from an initial phase space point Xq = 
(xo, 6*0) at the boundary between the system and the lead 
to the corresponding point X = (x, Q). The average of P 
over an ensemble or over energy gives a smooth function 
that reads 



(P(X,Xo;t)) 



cos(6') 



-t/TD 



(23) 



with the escape rate t^^ given in Eq. 

Using Eqs. {21]), and ([23]), we recover the Drude 
admittance 



3a/3 (^) 



where = 



M 



aOaP 



N 



1 



1 — iuJTY) 



(24) 



C. Weak localization for transmission, reflection 
and coherent backscattering 

1. Weak localization 

The leading-order weak-localization correction to the 
conductance was identified in Refs. [33[39| as those aris- 
ing from trajectories that are exponentially close almost 
everywhere except in the vicinity of an encounter. An 
example of such a trajectory pair for chaotic ballistic 
systems is shown in Fig. [21 At the encounter, separat- 
ing the 'loop' from the 'legs', one of the trajectories (7') 
intersects itself, while the other one (7) avoids the cross- 
ing. Thus, they travel along the loop they form in op- 
posite directions. In the semiclassical limit, only pairs 
of trajectories with a small crossing angle e contribute 
significantly to weak locaHzation. In this case, each tra- 
jectory remains correlated for some time on both sides 



of the encounter. In other words, the smallness of e re- 
quires two minimal times: Ti,{e) to form a loop, and 
Tw(e) in order for the legs to separate before escaping 
into different leads. The encounter introduces a typical 
length scale 5r±^ that corresponds to the perpendicular 
distance between the two paths in the vicinity of the 
encounter. In the case of hyperboHc dynamics, we get 
6r±^ = v-pe/(2X) ^ Le. Hence, the typical minimal time 
is given by T({e) = X'^ ln[{£/Sr±f], with £ = {L, W} 
that we can approximate as 

me) - X-Hn[e-% (25a) 
Tw(e) ~ X-Hn[e-'-{W/Lf]. (25b) 

The presence of the external driving does not change 
this picture. Each weak-localization contribution accu- 
mulates a phase difference given by the linearized action 
SS{Ef) ~ SSns = E^e^ />w^. Following the same lines 
as for the derivation of the Drude contribution, though 
the sum over paths is now restricted to paths with an 
encounter, the sum rule ((22|) still appHes, provided the 
probability P(X, Xq; t) is restricted to paths which cross 
themselves. To ensure this we write 

P(X,Xo;t) = / dR2dRiP(X,R2;t~t2) 

X P(R2,Ri;t2 -ii)P(Ri,Xo;ti), (26) 

where the integration is performed over the energy sur- 
face C. Here, we use R^ = (r^j^i), e [— 7r,7r] for phase 
space points inside the cavity, while X lies on the lead 
surface as before. 

We then restrict the probabilities inside the integral 
to trajectories which cross themselves at phase space 
positions Ri,2 with the first (or second) visit of the 
crossing occurring at time ti (or t2). We can write 
dR2 = w| sinedtidt2de and set R2 = (ri,0i ± e). Then 
the weak-localization correction is given by 



1 



-J dXoJ deSRe 









with. 



(P(Xo,e,c^)), 
(27) 



2z;|sine / dt dt2 / dh 

xpFCOs6io dY dRiP(X,R2;t-t2) 
Jr Je 

xP(R2, Ri; t2 - ti)P(Ri, Xo; h) e''"*. 



(28) 



Under our approximation 



tj = t, the intro- 



duction of the driving frequency leads to performing a 
Fourier transform of the survival probability, and we ob- 
tain 



(F(Xo,e,c^)) 



{vptd)'^Pf sin e cos 6q Na 



TTfl N 

exp [-rL/Tp] exp [w(Tl + Tw)] 



(29) 
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with n the cavity area. Inserting Eq. I|29|) into Eq. (|27|) . 
the e integral is dominated by small angle (e ^ 1) 
contributions, allowing for the approximation sine ~ e 
and pushing the upper limit to infinity. This yields 
an Euler Gamma function times an exponential term 

e-^E/rDg»^(rg'+r°P) (^j^j^ ^op ^cl gj^^j^ (fT|2|) 

that reads, to leading order in (Atd)^^, 



de2 5Re 



exp 



mVpTD 



Xh 



w 



Atd 



(30) 



Performing the Xq integral and using Np = {■7Th)^^pFWi3 
and TV = {hTu)^^mn, the weak-localization correction to 
the unscreened admittance is 



9ai3 y^j ~ 



(1 - 2za;rD) e'"(^E +^e'') 
(1 - iujTu)^ 



.(31) 



We note that due to the absence of unitarity of the un- 
screened admittance we need to explicitly evaluate all the 
elements of g'^p{uj). The weak-locaHzation contribution 

to reflection r^^a^iuj) is derived in the same manner as 
g^^'(w), replacing however the factor Np/N by Na/N. 
We then obtain 



-Tg'/rp (1 - 2ztJTD)e' 



(1 - iujTY)y 



(32) 

However as in the dc-case there is another leading- 
order contribution to the reflection, the so-called coherent 
backscattering. This differs from weak locaHzation as 
the path segments that hit the lead are correlated. This 
mechanism should be treated separately when computing 
the Ehrenfest time dependence, which is the object of the 
next paragraph. 



S. Coherent backscattering 

Though the correlation between two paths does not 
influence the treatment of the external frequency, it in- 
duces an action difference 6S{Ef) = SSchs = —{pn± + 
™Aroj^)roj_ where the perpendicular difference in po- 
sition and momentum are rQ± ~ {xq — x) cos 9o and 
Po± = — pf(^' — ^o)- As for weak localization, we can 



identify two timescales, ^T^,^T^, associated with the 
time for paths to spread to L, W, respectively. However 
unlike for weak localization we deflne these timescales 
as times measured from the lead rather than from the 
encounter. Thus we have 

2 

T^{rQ±,pQ±) — 1^ In [("T-^^)/ \Po± + mXroj_\] , (33) 

with £ — {L, W}^. Replacing the integral over Xq by 
an integral over (ro_L,po±) and using ppcos^odXo = 
dpQ±drQ±, the coherent-backscattering contribution 
reads 



(F=^^(Xo,c.)), 



J a 

with 

{F'^''%Xo,Lj)) ^ I dt /dXP(X,Xo;i)e 



(34) 



N 



1 — ILOTd 



(35) 



As in the dc-case we perform a change of variables 
Po± = Po± + mAro^. Then we push the poj. integral 
limit to inflnity and evaluate the roj_ integral over Wa- 
This result, 



dpo_ 



hs\n(pQi_Wa/h) 



Po± 



Po± 



(1-2.^xd) 



mXL 

TT/ie^^e'^^^E+^E'i + [(Atd)"^] , 
together with Eq. |[35| and Eq. ([34]) yields 



(1 - iujTi)) 



(36) 



(37) 



Surprisingly the coherent-backscattering contribution 
thus has exactly the same exponential dependence on 
and as the other weak- locaHzation contributions. 
While in the dc-case this property is a consequence of 
current conservation, this fact is not obvious in the ac- 
case. 

At this point we can summarize our results for the 
unscreened admittance. Prom Eqs. ^ [311 [32l \Tl\ . 

g^p{oj)) can be written as 



cxp 



CXp ICiJ^T] 



N{1 - iLUTn) 



N{1 



iLUTu) 

I 



A^ff(l-2zcjTD) 
iV(l - iajT^y 



-5^p \+0{N-'). (38) 



Pirst we note that in the limit of zero Ehrenfest time we recover the RMT result for the unscreened admit- 
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tance of Brouwer and Biittiker— . Concerning the Ehren- 
fest time dependence of the admittance, we note that the 
result is consistent with the absorption study performed 
in Ref. [g^I- As for the dc-case we find the absence of 
the Ehrenfest time r^^ in the term exp[— /td] which 
derives from the classical correlation between the paths 
that constitute the encounter. The physical origin of the 
term exp [^(rg + t^)] comes from the fact that both 
trajectories that contribute to weak localization and co- 
herent backscattering involve an encounter that has a 
minimal duration of (r^ + t^^) (Leg part and loop part 
of the encounter, see Fig. [2]). The presence of this mini- 
mal duration, 2Tg = + t^^ , is in accordance with the 
Ehrenfest time shift prediction of the quantum correction 
to the survival probability^ and the photofragmentation 
statistics^. We return to the Ehrenfest time dependence 
in Sect. HVEl 

We can also consider the effect of a magnetic flux on the 
mesoscopic admittance. A weak magnetic fleld has little 
effect on the classical dynamics but generates a phase 
difference between two trajectories that travel in opposite 
directions around a weak-localization generating closed 
loop. This phase difference is <&/$o, where $o is the flux 
quantum, and <i> is proportional to the flux through the 
directed area enclosed by the loop. To incorporate this in 



the previous semiclassical treatment we must introduce a 
factor cxp[i$/$o] into F in Eq. ^ and F"^"" in Eq. ((Ml)- 
The calculation gives a Lorentzian shap o^^i^^i^^ for the 
^-dependence of the quantum correction to the average 
admittance, 



ti.wl/cbs/ 



u.wl/ cbs / ri\ 
.9a/3 ('^.0) 



1 4- ^1^$^ (rf/TD - iuJTfY 



(39) 



Here = afi^, with a a system dependent parameter of 
order unity, fl the cavity area and T[ is the time of flight 
between two consecutive bounces at the cavity wall. 



D. The screened admittance 

Following the self-consistent approach, the screened 
admittance is straightforwardly obtained when we sub- 
stitute Eq. ((38|) into Eq. ifTS]) and expand the result to 
leading order in N~^. This simple substitution is justi- 
fled, because the typical fluctuations of the unscreened 
admittance are of order N~'^. The screened admittance 
then reads 



N{l-iiOT) 



exp 



cxp [iuj{T§^ + T°^)] 



N{1 - iwtd) 



iV^(l - 2iLUT) 
N{1 - iujtY 



-5c.p\+0{N-'), (40) 



where = r^^ +NGo/C is the charge relaxation time 
or quantum RC time. Eq. l(40|) is the first intermediate 
result from which we can draw some general conclusions. 
At zero Ehrenfest time we recover the two-terminal result 
of Brouwer and Biittiker in Ref. js^l ■ The comparison be- 
tween the screened (Eq. (|40l) ) and unscreened (Eq. l(38|) ) 
admittance shows that the screening amounts to the re- 
placement of the dwell time td by the RC time r every- 
where up to the prefactor of the third term. Only for the 
weak localization and the coherent-backscattering contri- 
butions does the dwell time dependence survive. Though 
the relevant time scale for the classical admittance is the 
charge relaxation time r, the quantum corrections are 
characterized by the dwell time td. It is important to 
remember that td is a characteristic time scale of the 
non-interacting system. Its relevance here has its origin 
in the fact that weak localization is due to the interfer- 
ence of electronic waves, which is unimportant for charge 
accumulation in the system. The absence of the RC time 
T at leading order in uj is thus quite natural. We recall 
that, as constructed in the framework of the model, the 
admittance matrix Eq. (fjO)) is current conserving if the 
gate is included. The elements of the admittance related 
to the gate are obtained via the sum rule ^ . Neverthe- 



less, if we impose this above sum rule to the unscreened 
result we also obtain a conserved current, and this situ- 
ation corresponds to a cavity which has infinite capaci- 
tance to the gate. In the reverse limit of zero capacitance 
we reach the charge neutral regime that corresponds to 
putting T = in Eq. (|40)) . Upon performing that, we 
recover the charge-neutral limit obtained by Aleiner and 
Larkin in Refs. [S^fCTj which for the conventional weak- 
localization contribution reads 



wl,r— 0/ \ 
9al3 (^) 



Ar2 



exp 



(1 - iljtd) 



(41) 



We note that for the particular geometry of a capacitor 
(only one lead and one gate), since Eq. l(40|) is valid for 
any capacitive coupling, we can obtain the effect of the 
Ehrenfest time scale on the interference correction to the 
admittance of a mesoscopic capacitor. This was not pos- 
sible within the charge-neutral limit approach of Aleiner 
and Larkin, since the interference corrections considered 
here are absent in that case. 

Here one important remark is due. In both, Eq. l(40|) 
and Eq. (j4l|) the admittance involves an oscillatory be- 
havior as a function of the Ehrenfest time, which should 
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in principle be more easily accessible experimentally. In- 
deed, we see here in our quest for the Ehrenfest time 
physics a clear advantage in investigating weak localiza- 
tion in the ac-regime. In the static case, the ratio te/td 
is the only relevant and tunable parameter for the dc 
weak-localization correction. Consequently, the range of 
experimental investigation is considerably reduced by the 
logarithmic dependence of te on the system size. For 
the dynamical weak localization the frequency depen- 
dence w combined with the capacitive coupling C pro- 
vides more freedom in probing TE-behavior. However, 
although the ujte Ehrenfest time dependence was pre- 
dicted in Ref. [s^l (in which some possible experimental 
verification was forecasted in a magnetocondutance ex- 
periment or in an optical backscattering experiment) , we 
are not aware of any experimental verification of the ex- 
istence of such an oscillation. To date there exist only 
two experiments devoted to exploring the te signature: 
The shot noise experiment by Oberholzer et al.— and 
the weak localization experiment in an antidot lattice by 
Yevtushenko et al.— . Both experiments were performed 
in the static case. 



Here we recall that 2r| = t^' . At zero Ehrenfest 

time, Te = 0, we see as in Ref. |57f that while the initial 
time dependence is determined by td (first term of rhs of 
Eq. l(43|) ). for times larger than t* the t^-term in the log 
will be important. We therefore find a deviation from 
the classical exponential behavior. 

This conclusion still holds at finite Ehrenfest time, up 
to the inclusion of a time shift 2rg as predicted in the 
recent semiclassical derivation^ of the survival probabil- 
ity. 

The treatment of the screened case is more demanding 
due to the presence of the RC time t. However since the 
pole linked to the dwell time td is only simple, it is clear 
that even at incomplete screening, there is no term pro- 
portional to i^. This is in accordance with the absence of 
deviations for the interacting admittance. However, the 
Ehrenfest time dependence will be equivalent to the un- 
screened one, leading to a time shift . Only for complete 
screening (t = 0) it is possible to obtain a simple result, 
which reads on a log scale 



In 



TD 



In 



1 

TV 



(44) 



E. Pulsed cavities 

In this section we comment on the Ehrenfest time de- 
pendence of the admittance and its link to that of the 
survival probability^"^. To this end we consider the par- 
ticular case of a pulsed cavity^, i.e. the application of a 
pulse Xiaii) = aa5(t) to one of the contacts a. The re- 
sponse current at contact (3 to such a pulse will be propor- 
tional to the frequency integral over the ac-conductance. 



1 

2^ 



3a/3('^)cxp(-iwt). 



(42) 



This problem was previously addressed in Ref. [57| where 
the connection between the RMT calculation of the ad- 
mittance and RMT results for the quantum and the clas- 
sical survival probability22i2i were discussed. More pre- 
cisely, in Refs. |70ll7l| a difference between the quantum 
and the classical survival probability was predicted for 
times of order t* — ^tdth. The conclusion of Ref. jH3| 
was two-fold: first, based on the weak-localization cor- 
rection, a deviation of the unscreened admittance at t* 
was confirmed, while secondly the screened system was 
shown not to exhibit such a t*-dependence. 

Based on our semiclassical results p8|40p we are able 
to confirm this dependence. For the unscreened admit- 
tance, the weak-localization and coherent-backscattering 

contribution, 5g'^f^{t) = gap\t)+9a(t''{^)-' yields a com- 
plicated time-dependence and reads on a log scale 



In 



iVTD 



t - T°P 



TD 



(43) 



+ ln 



1 

N 



t_2Tl 

TD 



2td 



V. MULTI-TERMINAL SYSTEM WITH 
TUNNEL BARRIER 

The calculation of the admittance with tunnel barriers 
follows the trajectory-based method recently developed 
by Whitney^ for the dc-case. We recall here the three 
main changes in the theory with respect to the trans- 
parent case. For more details on the inclusion of tunnel 
barriers we refer to Ref. (45| . 

At first, in the presence of tunnel barriers the com- 
plex amplitude in Eq. lfT7|) is extended to include the 
tunneling probabilities reading^. 



C7 tp.ita.j 



Tfi- 



l3^T (/?',/) 



(45) 



where = |(dpa;o/da;)-y| is the rate of change of the ini- 
tial momentum for the exit position x of 7, 3\f-y(/3', j') 
is the number of times that 7 is refiected back into the 
system from the tunnel barrier on lead /3' and the trans- 
mission and refection amplitudes at the lead (i satisfy 
1^/3, = (1 ~ k/3,iP) = r,3,i- We note that without any 
loss of generality, we associated in Eq. (|45| the momen- 
tum (or px) with the channel i (or j). 

At this point the replacement of the semiclassical am- 
plitudes by their corresponding classical probabilities still 
holds, though the tunneling probabilities are included. 
As an example the probability to go from a phase point 
Xo (here we associate the channel i to the momentum 
ppcos^o) on lead /3 to an arbitrary point on lead a sim- 
ply satisfies (for a ^ P), 



[ dt ( dX(P(X,Xo;t)) 

JO J a 



(1) 



(46) 
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Encounter see Fig. 2 
leg 

Figure 3: A failed coherent-backscattering contribution to ac- 
conductance, g^'^^"" (uj) . It involves paths which return close 
but anti-parallel to themselves at lead a, but are reflected off 
the tunnel-barrier, remaining in the cavity to finally escape 
via lead /3. The cross-hatched area denotes the region where 
the two solid paths are paired (within Wa — W of each other). 



where we let T 



(1) 



/3 - Ef=fi and define 3^ = . 
More importantly, the introduction of a tunnel barrier 
induces three changes: (i) The dwell time (single path 
survival time) becomes 



and define 3\f : 



'Dl 



(47) 



because a typical path may hit a lead but be refiected 
off the tunnel barrier (remaining in the cavity) numerous 
times before tunneling and escaping. 

(ii) The paired-paths survival time for paths closer 
than the lead width is no longer equal to the dwell time 
instead it is given by 



'D2 



(48) 



where v'i^ = and we define H = Ea Ta ^. This 

is because a second path following a path which has not 
escaped will hit the same tunnel barrier, and thus may 
escape even though the first path did not. Compare this 
with a system without tunnel barriers: there a path has 
not escaped because it has not touched the leads; thus a 
second path following the first one has no possibility to 
escape. 

(iii) The coherent backscattering peak contributes to 
transmission as well as refiection. The positive contri- 
bution to the transmission competes with the usual neg- 
ative weak-localization contribution to transmission, see 
also Figd 

For the calculation of the Drude conductance, only 
change (i) above is required, yielding 



i R Oaf3 



r(i)r(i) 

i- Q J- ^ 



1 



74 1 — ILUTDl 



(49) 



When calculating the conventional weak-localization 
contribution we need changes (i) and (ii) above. Since 
the classical paths considered stay close to itself for a 
time Tvi/(e)/2 on either side of the encounter we must 
use the paired-paths survival time, rD2, for these parts of 
the path. Elsewhere the escape time is given by the single 
path survival time, tdi. With these new ingredients we 
find that the conventional weak-localization contribution 
becomes 



9aP 



(1 - iujTui)^ 



(50) 

with 8t-e = t'^ /t-d2 + (■''e ^ Te')/'''^^- The exponential 
suppression exp(— Ore) related to the classical correla- 
tion is simply the probability that the path segments 
survive a time as a pair {t^ /2 on either side of 
the crossing) and survive an additional time (r^ — r^^) 
unpaired (to complete a loop of length ). Similarly 
as for the transparent case, the exponential dependence 
exp[za;(rg +Tg^)] indicates that the minimal duration of 
a weak localization trajectory is -|- t^^. 

However as realized by Whitney^, this is not the to- 
tal weak-localization contribution to conductance, be- 
cause of failed coherent-backscattering g^'^^^{uj) that 
contributes to conductance (change (iii) above). We re- 
call that this involves a path which returns close but 
anti-parallel to itself at lead a, but is then refiected off 
the tunnel-barrier on lead a, remaining in the cavity un- 
til it eventually escapes through lead An example of 
such a trajectory is shown in Fig. [31 We can calculate 
the backscattering contribution as before but using Tt)2, 
when the paths are within Wa of each other, and tdi 
elsewhere. This result is then multiplied by the proba- 
bility that the path reflects off lead a and then escapes 
through lead /? and weighted by the dynamical factor 
(1 — iujTui)^^ due to the diagonal transmission from a 
to (3 i.e. the leg part of Fig [3l In addition to the co- 
herent backscattering expression for r'^'^^^[Lo) this gives 
a contribution to the admittance of the form 



tt,cbsl/ \ 
5a/3 i^) = 



u,cbs2/ \ 



u.chs / \ 



p{2) pi^j 



(1 - icjrDi)2 



(1) p(i) 



,(51a) 



.(2) 



^ /3 i a 



(1 - iLUTui)^ W 



Sal3 



p(2; 



1 — iujT-Di 'N 



(51c) 



where we recall that Ta'^ = EilTi i- 

Using Eqs. I|491 [SO j ISTj) . the unscreened admittance in 
the presence of tunnel barriers reads 
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^ a 



.(1) 



p(i)p(i) 



(1 - iujTDl) 



(1 - icJTDl) 



+ 



p(2) 



(1 - iuJTBl) 




p(i) p(i) 



As a check of the formula (|52|) . we can easily recover 
the previous Eq. (|38|) for the unscreened admittance ob- 
tained for transparent barriers and also the tunnel dc- 



conductance^. 

After the substitution of Eq. ([52l) into Eq. (H!]) the 
screened admittance in presence of tunnel barriers reads 



p(i)p(i) 

(i)r(i) _,:„f^£>,,°pW, r(2) /-p(i) ^ -p(2) _ o 
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N2 (l-iwTDi) I {l-iujTY {1-iujt) rL^) 




where the quantum RC time reads now r^^ = t^I + 
'NGq/C. We emphasize that from Eq. l(53|) it is possi- 
ble to derive all the results presented in this paper and 
therefore this equation is the central result of this paper. 

In the second line of Eq. l(53|) . the second contribu- 
tion in the brackets represents the correction due to the 
presence of the failed coherent backscattering. Impor- 
tantly, Eq. I|53p includes both, the limit of infinite capac- 
itance C and the transparent case. In the charge neu- 
trality limit (r ~ 0) the presence of the tunnel barriers 
does not drastically alter the conclusion drawn for the 
transparent case. Indeed, for the weak-localization cor- 
rection, in addition to the expected substitution Na,N 
by ri'\7^ , we observe only a renormalisation by a factor 
(ri^Vr^ ^ + r^^'/r^^^ - ^/^). Thus Eq. gl]) becomes 

5S^=°(-) = (54) 

l^rW + r« 1^) DV2 ■ 

More importantly, one of the main effects of the tunnel 
barrier in the dc-case was the suppression of the weak- 
localization correctior4^^ for opaque barriers. This 
suppression results from the competition between two 
purely quantum effects, interference and tunneling. The 
corresponding semiclassical treatment^ shows that the 
cancellation is due to an exact compensation between 
the weak-localization correction and the failed coher- 
ent backscattering. It is interesting that this conclu- 
sion cannot be generalized to ac-transport. Since the 
frequency dependence of the weak-localization correction 
differs from the one of the failed coherent backscattering 



the compensation cannot occur. Dynamical weak local- 
ization is thus more robust against the presence of tunnel 
barriers. We note, however, that for r = we recover the 
cancellation of the weak-localization correction with tun- 
nel probabilities, see Eq. I|54p . 

VI. CHARGE RELAXATION RESISTANCE OF 
A MESOSCOPIC CHAOTIC CAPACITOR 

To illustrate and apply the general results derived 
above, we consider here the mesoscopic equivalent of a 
classical RC circuit^^. A quantum coherent capacitor has 
been recently investigated experimentally by Gabelli et 
al.— using a two-dimensional electron gas. The quantum 
capacitor is composed of a macroscopic metallic electrode 
on top of a lateral quantum dot defining the second elec- 
trode. The role of the resistance is played by a quan- 
tum point contact that connects the quantum dot to a 
reservoir. The experiment was performed in the coherent 
regime at high magnetic field in the one edge state limit. 
Measuring the real and imaginary part of the admittance 
of such a circuit, Ref. ^ confirmed the predicted^^ uni- 
versal value of the quantized charge relaxation resistance 
of a single channel cavity, which is equal to half a resis- 
tance quantum h/2e^. 

Based on this experimental realization we propose here 
to investigate the opposite regime of large channel num- 
bers at zero magnetic field. This regime is not character- 
ized by the universal value of the preceding fully quantum 
one, however it should be experimentally accessible. If 
we assume that the quantum dot is chaotic we can map 
this system to the one-terminal geometry of the more 
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general set-up considered in the previous section. The 
transparency of the quantum point contact is replaced 
by the transmission probability of the tunnel barrier Fi. 
To simplify the result we assume in the following that the 
N channels of the capacitor have the same tunnel rate, 
i.e. Fi i = r (Vi), the dwell time of the capacitor is thus 
TD = rH/(A^F). 

In a quantum coherent capacitor, there is obviously no 
dc-current, but we can address ac-transport via the ad- 
mittance G{ui)2^^. At low temperatures it is character- 
ized by an electrochemical capacitance and a charge 
relaxation resistance i?„ 



"-<?) 



-iujCn 



LO^ClRq 



(55) 



In contrast to their classical counterparts, and Rq 
strongly depend on the local density inside the sample^. 
They are thus sensitive to the phase coherent dynamics 
of the electrons inside the sample and thus subject to 
dephasing. 

Uo{lo),Iq{uj) 
Gate 



The unscreened admittance elements are given in 
Eq. (|52| . The survival times tdi and td2 of this two-lead 
geometry are related to the real dwell time td of our 
capacitor and to the dephasing time = th/(A^0F0), 
where N^f, and F^^^ = F^ (Vi) are, respectively, the num- 
ber of channels and the tunneling rates of the dephasing 
leadi^: 



Tdi = Td 



td2 = Td 



1 + ^ 

^0 



(2-F) 



TD 

T0 



2-Fw 



(58a) 
(58b) 



Inserting expression ((52|) for the unscreened admittance 
elements into Eqs. l(56 t [F7 |) and performing an expansion 
in uj we get 



GoRq 



1 , i?(r^',r°P,r^) 



FiV 



(F7V)2 



(59a) 



Lead 1 



1 

1 



r 



Dephasing Lead 



Figure 4: Schematic picture of the mesoscopic capacitor with 
the dephasing lead (see text). The chaotic cavity has an extra 
lead (lead 0), whose voltage is chosen to render the net current 
zero, which leads to dephasing without a loss of particles. 
Since each channel has the same tunnel rate Fi.i = F and 
F^,! — r^, Vz, the dwell time of the capacitor is td oc {Nr)~^ 
and the dephasing time reads oc (A^^F^)"^. 



To model the loss of coherence of electrons inside the 
cavity we appeal to the so-called voltage/dephasing probe 
model^i, which consists of adding another lead (j), (see 
Fig. H]) to our cavity and tuning the potential of this 
probe in such a way that the net current is zero. Conse- 
quently any electron that enters this lead is immediately 
replaced by another one with an unrelated phase lead- 
ing to incoherence without loss of particles. While such 
an approach has recently been used for the mesoscopic 
capacitor in the one channel Hmitji^, here we investigate 
similar effects of the dephasing in the limit of large chan- 
nel numbers^^ where our semiclassical method is fully 
justified. 

The admittance can be written as 



G(c.) = 



-iuoC + x{^) ' 



where 



X{uj) = Go 9^1 ('^) 



(56) 



(57) 



where we additionally used the relation between the mean 
density of states, v, and the dwell time, td = hv/ {dsNT). 
The dephasing function D{t§^,t^^ ,t^) reads 



op ci 
i^(r-,.-,r,).Fe-^(--)-^&^ 



1 + ^ 

(60) 

We finally consider the effect of a magnetic fiux on 
the charge relaxation resistance. Substituting Eq. l[39|) 
(the dwell time being replaced by the survival time tdi) 
into Eq. l(56l) leaves the electrochemical capacitance 



unchanged; only the dephasing function D{t§^, t^, t^) is 
affected and replaced by 



(i-r) 



^(T^',T°P,r^,$) = Fe ^5 ^ 



1 + yi2$2iD ID 



(61) 

From this semiclassical investigation of the charge re- 
laxation resistance, we can see that the fully coherent 
limit (r0 = oo, = = 0, $ = 0) delivers. 



11/. 1 
""^ ~ G^TN \ ' ■ N 



Ra = — — (1 + ^1 +0(iV-3). 



(62) 



Eq. (|62| is the first derivation of the charge relaxation 
resitance in the large limit in presence of tunnel bar- 
riers. While the leading order was guessed^, the weak- 
localization correction to Rq has never been calculated 
before. Surprisingly, it is linear in the inverse tunnel rate 
r~^, indicating that the calculation of the sub-leading or- 
der correction cannot be simply obtained by an effective 
renormalisation of the channel number N^g = TN. 
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For the incoherent Hmit, obtained either by te — > oo, 
$ — > cxD or T0 = 0, we get a suppression of the weak- 
locaUzation correction and thus Rq reduces to 



This value corresponds to the fully incoherent limit 
that corresponds to the two-terminal resistance, and has 
been obtained under the simple application of our de- 
phasing process. Interestingly, this limit was not trivial 
to obtain in the edge state calculatioi>i^ {N = 1), where 
perfect inter-channel relaxation inside the voltage probe 
was assumed. This seems not to be required in the fully 
chaotic case in the limit N ^ 1. 



VII. CONCLUSIONS 

In this work we focused on the topic of ac-transport 
through chaotic balhstic conductors, addressing in par- 
ticular weak locahzation corrections to the admittance 
from a semiclassical perspective. Employing trajectory- 
based methods we confirmed RMT results for the bare 
and screened admittance and, going beyond RMT, de- 
rived the Ehrenfest time dependence. The Ehrenfest 
timescale enters twice into the expressions for dynami- 
cal weak localization: first, as an exponential suppres- 
sion with an exponent given by the ratio of the Ehren- 
fest and dwell time, te/td; second the dynamical weak 
localization acquires an oscillatory frequency-dependent 
behavior of period 2te, which may be amenable to mea- 
surements based on variations of the ac-frequency. We 
emphasize that our results are valid for any finite capac- 
itance C and hence not limited to the electroneutrality 
assumption of Ref. [33|. This extends the class of exper- 
imental settings for which the Ehrenfest time correction 
can be investigated. More generally, the results presented 
underline, firstly, the power of semiclassical techniques 
to provide a clear and quantitative picture of ac-driven 



quantum transport in the various regimes and, secondly, 
they give a justification of the "stub model"— in the low- 
frequency regime. 

Moreover we took into account tunnel barriers in the 
semiclassical approach to the ac-admittance, extending 
the work of Whitney^ on dc-transport. This led us to a 
general formulation of ac-transport. One main conclusion 
is that weak-localization is more robust against effects of 
tunnel barriers in the dynamical than in the dc-regime. 
The extension of our semiclassical treatment to tunnel 
barriers also enables us to access the experimentally rel- 
evant case of a quantum coherent capacitor, for which 
we provide the first derivation of the weak-localization 
correction to the charge relaxation resistance in presence 
of tunnel barriers. 

We add that, from a methodological point of view, 
the semiclassical approach presented might be helpful to 
achieve a better understanding of the proximity effect 
on the density of states of chaotic Andreev billiards. Fi- 
nally, the ac-conductance discussed here is closely related 
to problems of computing (photo-)absorption and, more 
generahy, hnear-response based dynamical susceptibili- 
ties for mesoscopic quantum systems. It appears promis- 
ing to apply the semiclassical techniques, developed here 
for (ac-)quantum transport, to refine earlier semiclassi- 
cal approaches^ to (photo-) absorption in closed ballistic 
cavities or metal clusters, which additionally poses the 
chahenge to semiclassically cope with screening effects 
and plasmon excitations. 
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